Vortex Structure in Underdoped Cuprates 
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In underdoped cuprates the normal state is highly anomalous and is characterized as a pseudogap 
phase. The question of how to describe the "normal" core of a superconducting vortex is an out- 
standing open problem. We show that the SU(2) formulation of the t-J model provides a description 
of the normal state as well as the vortex core. Interestingly, the pseudogap persists inside the core. 
We also found that it is likely that the core consists of a state which breaks translational symmetry 
due to the existence of staggered current which generate staggered magnetic field with very slow 
dynamics. This staggered flux state is likely to be the ground state for magnetic fields higher than 
Hc2- Experiments to test this picture are proposed. 



I. INTRODUCTION 



It is now widely appreciated that high Tc supercon- 
ductors are fundamentally different from conventional su- 
perconductors in that they emerge by introducing doped 
holes into a Mott insulator. This contrast is most ap- 
parent in the underdoped region where the density x of 
doped holes is small. Experimentally, this is also the 
regime where the physical properties are most anoma- 
lous. 0] Much attention has been focused on the nor- 
mal state, which is characterized by a pseudogap regime 
below a relatively high temperature T* « iOOK . The 
pseudogap appears in spin excitations and in tunneling 
and ARPES experiments. The superconducting state is 
anomalous as well, in that the superfluid density is pro- 
portional to the hole density x and not the electron den- 
sity (Fermi surface area) 1 — a; as in conventional su- 
perconductors. Recently, it has become possible to per- 
form STM tunneling in the superconducting state and 
probe the electronic structure of the vortex core. 
This raises the following interesting question. Common 
sense would indicate that the vortex core should be made 
up of the normal state and one would expect the pseu- 
dogap, i.e. a dip in the tunneling density of states, to 
persist in the core region. This is in fact what is seen 
experimentally. Yet a conventional description of a vor- 
tex core requires that the order parameter vanishes inside 
the core, which is usually accompanied by the vanishing 
of the energy gap. Thus it is clear that the electronic 
structure of the vortex core in the underdoped region 
is qualitatively different from that given by conventional 
theory. This point was made eloquently in a recent paper 
by Franz and Tesanovic (FT). ||| 

It is clear that any attempt to model the underdoped 
vortex core must include the physics of the proximity to 
the Mott insulator, i.e. the strong correlation physics. 
One of the few analytic tools available for this purpose 
is the slave boson method used to treat the constraint of 
no double occupation in a strong correlation model such 



as the t-J model. FT employed the U(l) formulation of 
this theory, where the electron operator Cia is written as 
Cia = fiubl and the no double occupation constraint is 
replaced by fj^fia + blbi = 1, which is in turn accom- 
plished by the introduction of a U(l) gauge field a. In 
mean field theory the pseudogap gap state is described 
by a pairing of the fermions, A/(j7) = {f^ ,ifi,i+fj) , where 
r) is a nearest-neighbor vector and A(?7) has d symme- 
try- The superconducting state is described by bose 
condensation of the bosons (6) ^0. A/ is not gauge 
invariant and the onset of the pseudogap is merely a 
cross-over, but the appearance of (6) = bo ^ triggers 
the appearance of the superconducting pairing amplitude 
(c^iCii+fi) = b^Af which is gauge invariant and physi- 
cal. Within this theory, FT propose a description of the 
vortex state where the bosonic amplitude (b) vanishes 
inside the core but the fermion pairing amplitude |A/| 
remains finite. Since the electronic spectrum is given by 
the fermion dispersion, the core will retain the energy 
gap, just as in the pseudogap state. 

Upon closer examination, FT pointed out that this so- 
lution requires that the gauge field has negligible restor- 
ing force, i.e. a "Maxwell" term of the form cr(V x a)^ 
must have very small coefficient a. This requirement is in 
fact related to a problem discussed by Sachdev Q and by 
Nagaosa and Lee [|| some time ago. Due to the existence 
of two fields Af and (6), it is possible to construct sev- 
eral kinds of vortices. The field Ay is minimally coupled 

to a in the form | + ^a) A/|^, whereas the field (6) is 
coupled to a combination of a and the electromagnetic 
field A in the form \ {^ + a - ^A) (&)|^ 

The different kinds of vortices are as follows: 

(i.) A vortex carrying the conventional hc/2e fiux quan- 
tum. A gauge vortex carrying half a flux quantum 
i/i is generated so that (6) has no singularity. The 
phase of A^ winds by 27r and its amplitude vanishes 
in the core. This is just like the conventional vor- 
tex in that the energy gap vanishes inside the core. 
This describes the optimal or overdoped region. 
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(ii.) An hc/e vortex. This involves no winding of Ay 
and no gauge flux. The advantage is that |A/| is fi- 
nite in the core and the pseudogap is preserved. 
This state is energetically favorable because the 
cost of the boson vortex is small for small x. The 
price one pays is that because the boson carries 
charge e, this vortex carries a double superconduct- 
ing flux quantum This has so far not been ob- 
served. 

(iii.) The FT vortex. A third possibility proposed by FT 
is that a flux tube for the gauge field a carrying 
— \h gauge flux is attached to the core. Now the 
A flux can be a conventional flux quantum ^ and 
the phase of (&) winds by 27r, with (6) = in the 
core. On the other hand, A/ sees only a flux tube 
and remains non-zero in the core. It is this latter 
requirement which forces the a flux to be a flux 
tube, i.e. confined to a lattice plaquette. 

Actually this possibility was considered by Nagaosa 
and Lee and dismissed because the energy cost of a flux 
tube is large in the presence of a Maxwell term. The point 
is that the theory for A/ and (6) must be considered 
as a low energy effective action, and terms allowed by 
symmetry such as the Maxwell term will be generated 
by eliminating the high energy degrees of freedom. We 
expect the energy of the flux tube will be of the order of 
the cut-off scale, i.e. the fermion band width J. This will 
make this kind of vortex very costly in energy compared 
with the hc/e vortex in the limit of small x. 

FT appealed to the papers by Nayak ||^ and D.H. Lee 
1^ to justify setting cr = 0. Even assuming for the sake 
of argument that a vanishes and the a flux tube costs no 
energy, the FT vortex still had a core energy at least of 
order J. Although the pairing fleld A/ cannot see the 
h/2 flux tube of a, the fermions see the flux tube. The 
mismatch of phase tt at the lattice scale in the fermion 
wave function will cost an energy of order J . (Actually, 
this is why the h/2 flux tube costs an energy of order J 
as discussed above.) In order to reduce this energy cost, 
A f likes to vanish in a region of size of coherent length 
S^F ~ 'i'//^o, where Aq is the spin gap. Such a vortex 
has a core very similar to the standard BCS vortex (the 
case (i) mentioned above). The fermion contribution to 
the core energy is reduced to a value of order Ao. 

We should add that recently Senthil and Fisher ||l^] 
proposed a model of the vortex based on their Z(2) gauge 
theory which carries ^ flux quantum and contains a 
pseudogap in the core. This is accomplished by attaching 
a Z(2) vortex to the core. Senthil and Fisher |lj] recently 
showed how the Z(2) gauge theory can be placed in the 
context of the U(l) theory and it becomes clear that their 
model of the vortex is intimately related to that of FT. 
Senthil and Fisher combine the phase of the boson and 
half that of A/ to form the phase of the "chargon" which 



bose condenses. The Z(2) vortex is then the residue of the 
half flux tube of FT. The Z(2) vortex is also localized to 
a lattice plaquette and has an energy gap which Senthil 
and Fisher identify with the pseudogap scale. This also 
renders this vortex costly in comparison with the hc/e 
vortex (where the chargon winds by 47r) in the limit of 
small X. Thus we conclude that models based on U(l) 
mean-field theory still have difficulties coming up with a 
stable hc/2e vortex with a pseudogap core in the limit of 
small doping. 

Several years ago, we introduced an alternate formula- 
tion of the constraint in the t-J model called the SU(2) 
theory. This model is designed to connect smoothly 
to the Mott insulator at half-filling, in that the SU(2) 
symmetry known to be present at half-filling is preserved 
for finite doping. The SU(2) mean-field theory should 
have a better chance of describing the small doping limit. 
In this paper we show that this theory leads naturally to 
a stable hc/2e vortex in the underdoped limit. The spin 
gap is finite both inside and outside the vortex core. Pos- 
sible experimental consequences are explored at the end 
of the paper. 

II. REVIEW OF THE SU(2) FORMULATION 

First we summarize some of the salient features of the 
SU(2) formulation. This is well understood in the un- 

doped case, where SU(2) doublets Vtj = (/Tj'/|j) 

V'ii = ijiii '^6^6 introduced on each site j to rep- 

resent the destruction of spin up and spin down in the 
subspace of one fermion per site. Wen and Lee ex- 

tended the SU(2) formulation away from half filling by in- 
troducing a doublet of bosons hj — [bij , b2j ) . The phys- 
ical electron is represented as an SU(2) singlet formed 
out of the fermion and boson doublets: c^j — -^h^jip^j. 
The constraint of no double occupation is enforced by 
projecting onto the SU(2)-singlet subspace of the ex- 
tended hi,tpai Hilbert space. On each site there are 
three such singlets, corresponding to |spinup) = /||0), 

I spin down) = /||0) and 




The role of the two bosons can be visualized as follows. 
In contrast to the U(l) formulation, the fermions may re- 
main at half-filling upon doping. Then a typical fermion 
configuration will contain spin-up or spin-down singly oc- 
cupied sites, as well as empty and doubly occupied sites. 
The latter sites are both spin singlets and have the cor- 
rect spin quantum number for a vacancy. The hi boson 
is used to mark the empty site and the 62 boson the dou- 
bly occupied site, and both bi and 62 carry unit charge. 
This picture is a bit over simplified, in that it is a linear 
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superposition given by Eq. (1) which correctly specifies 
a physical hole. 

In order to perform the projection to SU(2) singlet, 
three sets of gauge fields a^j, associated with the three 
Pauli matrices ,£ = 1,2,3, are needed. These are the 
generalization of the time component of the gauge field 
floj in the U(l) formulation. The exchange and hopping 
terms are decoupled to give the mean-field Hamiltonian, 



H = 



T Ipaj + 



j 



3 '"3 



<]k> 



u. 



The matrix Ujk is given by 



(2) 
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FIG. 1. Pictorial representation of the staggered flux state. 
The hopping integrals are complex in such a way that the sum 
of the phase angle around a plaquette yields a net flux which 
alternates in sign from plaquette to plaquette. This gives rise 
to circulating fermion currents on the bonds as indicated by 
the arrow. In the presence of hole doping and condensation 
of the bosons, circulating physical hole currents appear. We 
refer to this state as the staggered flux state. 
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and is enforced 



The hole density is {b\bi + 62^2) = 
by the chemical potential /i. The Lagrangian associated 
with Eq. (2) is invariant under the local SU(2) gauge 
transformation 



V'aj g]i^aj , hj g]hj , Ujk g]Ujkgk , 



(5) 



where gj = exp(iAj • t) is a space and r dependent 2x2 
matrix that represents an SU(2) group element. 

In Ref. ||l5| the SU(2) mean-field theory was worked 
out by making the approximation that a^j is indepen- 
dent of space and r. Of special interest is the pseudogap 
phase which occupies the low doping part of the phase 
diagram. (Note that the pseudogap phase is called stag- 
gered flux or s-flux phase in the SU (2) theory of Refs. fl^ j 
and [0. Despite its name, the s-flux phase in the SU{2) 
theory is translation invariant and has no staggered phys- 
ical magnetic field. In this paper, we will reserve the 
name "staggered flux phase" for the staggered flux phase 
in the U{1) theory, which does have staggered physical 
magnetic field. (l9) We use "spin-gap phase" to refer to 
what we previously called the staggered flux phase in the 
SU{2) theory.) In the spin-gap phase Cq ■ — with finite 



nearest-neighbor Xjk = X ^'iid Ajj. = ztA-^ for x or y 
bonds. This resembles the fermion pairing phase in the 
U{1) theory. However, due to the gauge symmetry given 
by Eq. (5), the same mean- field state can be constructed 
with the choice 



(6) 



This resembles the staggered flux phase in U(l) mean- 
field theory ||l^ because the hopping matrix elements 
are complex and the sum of the phase angle around a 
plaquette gives a flux which alternates in sign, as shown 
in Fig. 1. Equation (6) and the fermion pairing state Eq. 
(3) both give the fermion dispersion. 



1/2 



where 



e{k) = —2Jx (sin k^a + sin kyO) 
r]{k) = 2JA/ (sinkxa — sinkya) . 



(7) 



(8) 
(9) 



Due to our gauge choice, this dispersion is shifted by 
(5' 5) compared with the more conventional parameter- 
ization which has a maximum Aq = JAj at (0, tt), (tt, 0) 
and nodes at (±-|,±^). The boson dispersion is the 
same except that J is replaced by t. 

The breaking of translation symmetry shown in Eq. 
(6) is a gauge artifact, because Eq. (6) is mapped onto 
Eq. (3) which is translationally invariant by a site de- 
pendent SU(2) transformation. (l|] Thus in the SU(2) 
mean-field theory, the spin-gap phase includes fluctua- 
tions between the pairing state, the staggered flux state, 
and many other states in the U(l) formulation. Note that 
the mean-field ansatz given in Eq. (6) is itself invariant 
under a rotation. Thus the SU(2) symmetry has been 
broken down to U(l) in the spin-gap phase. 

In Ref. [|l8| this point of view was clarified and the 
approximation of constant ap improved by introducing a 
nonlinear a model description in terms of a slowly vary- 
ing boson field. The idea is that at low temperatures the 
bosons are nearly bose condensed to the bottom of the 
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staggered flux 




staggered flux 



FIG. 2. The isospin quantization axis I represents differ- 
ent states depending on its orientation. In the north and 
south poles, it represents the staggered flux states. These 
are two degenerate states with the current pattern shifted by 
one lattice constant. In the equator it represents the d-wave 
superconductor. Vectors connected by rotation around the 
z axis are gauge equivalent and represent the same physical 
state. 



boson bands and are slowly varying in space and time. 
On the other hand, the fermions have a short coherence 
length — Vf/Ao which is the lattice scale because 
Aq J/3. Then the fermions follow the local boson 
field and can be integrated out, after choosing an elq field 
which minimizes the action locally. The result is an effec- 
tive Lagrangian which depends only on the local boson 
field. It is convenient to choose the fermion mean field in 
the staggered flux representation given by Eq. (6), be- 
cause the symmetry breaking pattern from SU(2) to U(I) 
is manifest. The bottom of the boson band is at (f, f) 
and we write hj = hj exp [—i{jx + jy)^) and consider hj 
to be slowly varying. At low temperatures, lijhj = x and 
we write 



the vector /, as shown in Fig. 2. I pointing along z 
corresponds to the staggered flux state in the U(l) for- 
mulation. The polarization of the boson fleld generates 
a nonzero Oq, which corresponds to a shift in the chemi- 
cal potential for the fermions. This in turn converts the 
nodes at (f , f ) into a small Fermi surface pocket. This 
state breaks translational symmetry and is characterized 
by a staggered pattern of physical hole current distribu- 
tion as shown in Fig. 1. We shall refer to this state as 
the staggered flux state. / pointing along —z describes 
the same physical state except that the current pattern 
in shifted by a unit cell. On the other hand, I in the 
x-y plane corresponds to a d-wave superconductor state 
(with a finite chemical potential) which does not break 
translational symmetry. This is not obvious in the rep- 
resentation given in Eq. (6) and it is best seen by per- 
forming a space dependent SU(2) rotation according to 
Eq. (5) to the representation given in Eq. (3). Note that 
I pointing anywhere in the x-y plane is gauge equivalent 
and corresponds to the same physical d-wave supercon- 
ducting state. This is a consequence of the residual U(l) 
symmetry of Eq. (6) and any vectors I related by a 
(j) rotation represents states that are gauge equivalent. 
On the other hand, I pointing somewhere in between 
the north pole and the x-y plane are distinct states with 
both superconducting pairing and translational symme- 
try breaking. This state of affairs is summarized by an 
effective Lagrangian as derived in Ref. . It takes the 
form of an anisotropic 0(4) ct- model (jzip + |z2p = 1) 
coupled to gauge fields. For the purpose of this paper we 
restrict our attention to time independent variation and 
the Lagrangian takes the simplified form 



xt\D.iz\ 



x'^J 



-\ziz2\' + - {\zi\' ^ \z,\'y 

Cl C3 



^ (3)rr (3) 



where 



Zj2 



(10) 



(14) 



(15) 



where \zja\ = 1 and are parametrized by 



Zj2 = e'"e*2 sin - 
2 



(11) 
(12) 



The phase a is the overall U(l) phase which couples to 
the electromagnetic field. The angles and cj) are best 
visualized by introducing the isospin quantization axis I. 

I = z^TapZfj = (sin 6 cos 0, sin 9 sin (f>, cos 9) (13) 

i.e. 9 and </> are the polar angles of the quantization 
axis. The physical state depends on the orientation of 



is the covariant spatial derivative {j = x,y), is the 
spatial component of the a'^' gauge field, and ci and 
C3 are numerical constants of order 1. Since the SU(2) 
symmetry has been broken down to U(l) by Eq. (6), 
the a*^^) and a'^'^^ gauge fields are massive by the Higgs 
mechanism and have been ignored. The first term in Eq. 
(14) is the boson kinetic energy minimally coupled to 
the remaining a^^^ gauge field and the electromagnetic 
field A. The second term is a phenomenological term 
introduced to describe the difference in energy between 
the superconducting state and the staggered flux state, so 
that the quantization axis prefers to lie in the x-y plane. 
The third term comes from integrating out the fermion 
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degrees of freedom where Hjk is the fermion polarization 
bubble. In momentum space it is given by [[L2[ 



(16) 



i.e. it does not take the Maxwell form which would have 
been proportional to q^. 

At higher temperatures the anisotropy term (second 
term) in Eq. (14) is unimportant and the quantiza- 
tion axis is disordered. This corresponds to the spin-gap 
phase. At low temperature, the quantization axis picks 
out a direction in the x-y plane and at the same time the 
U(l) symmetry corresponding to a^^^ is spontaneously 
broken. This corresponds to the d-wave superconductor. 

Recently, quasi long-range correlations in the staggered 
current have been found in the Gutzwiller projected d- 
wave BCS wave-function |20[ and in the exact ground 
state of small samples. Such current fluctuations are 
very natural in the SU(2) theory, and are a consequence 
of fluctuations of the quantization axis / towards the 
north and south poles. We have suggested that these 
staggered current correlations may characterize the pseu- 
dogap state, but experimental detection of such fluctu- 
ating currents seems to be very difficult. Now we are 
ready to use this picture to describe the vortex in the 
superconducting state, and show that the staggered cur- 
rent fluctuations may slow down inside the vortex core, 
making its detection more hopeful. 



III. MODEL OF THE VORTEX CORE 

Our model of the vortex is the following. Far away 
from the core — I62I, but 62 = \/x zi changes its phase 
[a -\- <^/2) by 27r as we go around the vortex, while h\ = 
yjx z\ does not change its phase. The vortex contains ^ 
flux for the A field and \ flux for the a^'^^ field. From Eq. 
(15), 62 sees the sum of a^'^) and A, i.e. a unit total fiux, 
while h\ sees no net flux, so the winding we suggested is 
consistent. Note that the average phase a [see Eq. (11) 
and (12)] has a winding of tt, as appropriate for an ^ 
vortex. 

As we approach the vortex core, the amplitude of \zi\ 
must vanish to avoid a divergent kinetic energy from the 
first term in Eq. (14). Thus the center of the vortex core 

is represented by {h = -^2;, 0^ and is just the staggered 
flux state. As shown in Fig. 3, the quantization axis I 
provides a nice way to visualize this structure. It points 
to the north pole in the center of the vortex and lies in 
the equator far away, but its azimuthal angle winds by 27r 
as we go around the vortex. This is sometimes referred to 
as the "Meron" configuration, or half of a hedgehog. It is 
important to recall that I parameterizes only the internal 
gauge degrees of freedom d and 0, and the winding of 
by 27r has nothing to do with the winding of the overall 
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FIG. 3. Structure of the superconducting vortex. Top: 61 
is constant while 62 vanishes at the center and its phase winds 
by 27r. Bottom: The isospin quantization axis points to the 
north pole at the center and rotates towards the equatorial 
plane as one moves out radially. The pattern is rotationally 
symmetric around the z axis. 



phase a by TT around the vortex. To visualize the winding 
of both a and 0, it is necessary to go back to the (61, 62) 
representation. 

We can make a rough estimate of the vortex energy. 
Assume that the core size (size of the Meron) is ic and 
the size ofthea(3) flux is £a- There are four contributions 
to the energy. The flrst is the energy difference between 
the superconducting state and the staggered flux state. 
The main energy cost comes from the Fermi pockets. 
Assuming the area of the pockets to be x, we estimate 
an energy cost of P^x^/'^y/ J On the other hand, the 
Meron size cannot be smaller than 1/%/a; without costing 
too much kinetic energy. (In fact, the effective action 
is valid only for momenta q < ^/x since we kept only 
the first quadratic term.) The second term comes from 
the electronic supercurrent and is of order xtlii(X/£c), 
where A is the London penetration depth. The third 
term comes from the supercurrent associated with the 
a*^'^) gauge field, which is of order xt\n{£a/f^c), assum- 
ing £a > tc- Finally, the fourth contribution is from the 
gauge field action, the last term in Eq. (14). Setting 
q — in Eq. (16), we estimate this contribution to be 
^2|a(3)|VJAo/4 « VJ^o/L- The important point is 
that unlike the U(l) case, the gauge field is not confined 
to a flux tube, but can spread over a distance £a- We note 
that the supercurrent contributions depend logarithmi- 
cally on £a and £c, so that the main dependence comes 
from the flrst and fourth contributions. The staggered 
flux core size £c would like to be as small as possible, 
while the size of the gauge flux £a would like to be large. 
However, our estimate of the gauge flux energy should be 
cut off for q < X, because bosonic contributions will enter 
Eq. (16). Thus we conclude that the staggered flux core 
occupies a radius of x~^^'^ while the gauge field occupies 
a radius of x^^. The above estimate is very crude. The 
main purpose is to show that a standard hc/2e vortex is 
possible with a staggered flux core which does not cost 
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too much energy as a; — > 0. 

If we include effects of fluctuations, the size of the 
staggered flux core will very likely be bigger than the 
above estimate. One way to include the fluctuation ef- 
fects is through the following consideration. We have 
shown that due to the excitation of quasi-particle, the 
superfluid density is reduced in the vicinity of the vor- 
tex core. Q We have also shown that the quasi-particles 
carry current evp after including the fluctuation effects. 
p3[ In this case the superfluid density vanishes inside a 
radius of , which we identify as the vortex core. ||2^] 
This argument gives a lower bound on the vortex radius, 
which matches the radius la- Inside this radius the su- 
perconducting state loses phase coherence and becomes 
more costly in energy. Thus our earlier estimate may 
have over estimated the energy difference between the 
staggered flux state and the superconducting state inside 
the core and the staggered flux state may expand to oc- 
cupies the entire core of radius where the superfluid 
density vanishes. The important point is the topological 
structure of the vortex, which should be robust, while 
the details of the structure may be model dependent. 

One important consequence of the topological struc- 
ture is that there are two kinds of vortices, because the 
isospin quantization I can also rotate to the south pole 
at the vortex core. This just expresses the fact that the 
staggered flux state is doubly degenerate, with the stag- 
gered flux shifted by one unit cell. In the normal state 
these degenerate states fluctuate between each other, be- 
ing smoothly connected via the superconducting state. 
Inside the vortex core of the superconducting state, the 
topological structure of the vortex forbids such smooth 
fluctuations, and freezes in the staggered current pat- 
tern. Thus the vortex core is closely related to, but not 
identical to, the pseudogap state. 

Since the degrees of freedom in a vortex core is finite, 
the two possible staggered flux states inside the core can 
tunnel into each other. If the staggered flux core is as 
small as x"^^^, the tunneling rate can be as large as the 
spin gap. However, if the staggered flux core has a size 
of order x~^ (which is more likely), the tunneling will be 
reduced exponentially. Dissipation due to quasiparticles 
may further suppress the tunneling rate. Indeed, this 
problem is analogous to the tunneling between degener- 
ate two level systems coupled to a Fermi sea. There due 
to the orthogonality catastrophe the tunneling rate can 
scale to zero and the state completely frozen for strong 
enough coupling. In general such exponential tunneling 
rate is difficult to calculate, but we are hopeful that the 
dynamics will slow down sufficiently for the staggered 
currents to be measured experimentally. 

As the magnetic field is increased, the vortex cores 
eventually overlap at H = Hc2- The staggered current 
states overlap and it is reasonable to believe that the 
ground state should be a long range ordered staggered 
flux state, especially if the staggered flux core has a size of 



order The unit cell is doubled and the ground state 
is a Fermi liquid, with small Fermi pockets with area x. 
We predicted that Hc2 ^ x'^, since the core size scales 
as x"^. If a high quality underdoped sample can 
be made, Hc2 can be at a scale amenable to laboratory 
experiment. The Fermi pocket may be measurable by 
cyclotron resonance or Shubnikov-de Haas experiments. 
The cyclotron resonance has a unique signature because 
the Fermi surface is close to a Dirac point so that the 
Landau levels are not uniformly spaced. The doubling 
of the unit cell is difficult to measure directly, because 
the staggered current pattern does not couple to charge 
density modulation. It does produce a small staggered 
magnetic field, which we estimate very crudely to be of 
order 10 gauss. ||l9|,0 The possibility of detecting the 
staggered magnetic field by neutron scattering and fi- 
SR was investigated theoretically by Hsu et al. ||l^ They 
estimated the neutron scattering intensity to be 1/70 of 
that from the ordered moments in the insulator. 

IV. EXPERIMENTAL PROBE OF THE 
STAGGERED CURRENT 

If it is not possible to reach H > Hc2, the topologi- 
cal aspect of the vortex offers us an opportunity to test 
the staggered current picture. It is difficult to probe the 
staggered current pattern in the normal state because of 
spatial and temporal fluctuations. One of the few possi- 
ble techniques is X-ray scattering which couples to chi- 
rality fiuctuations at (7r,7r). However, according to 
our analysis, the dynamics of the staggered current pat- 
tern slow down inside the vortex core. Depending on the 
time scale, it may be possible to measure the small stag- 
gered magnetic field created by the circulating current. 
The field distribution in the vortex state is remarkably 
uniform, as expected for an extreme type II superconduc- 
tor. From /i-SR measurements, the field distribution has 
a width of roughly 5 gauss at H = 0.5T. |^ It should be 
even narrower at higher fields. If the dynamics is slower 
than the /x-SR scale, the field distribution inside the vor- 
tex core is detectable. For even slower dynamics, a more 
sensitive experiment is NMR. In YBCO, the Y ion is ide- 
ally placed to detect this current, because it sits at the 
center of the plaquette. The weak magnetic field gener- 
ated by the circulating current will produce side bands 
in the Y NMR line, with a splitting independent of H 
but with weight proportional to H. However, there re- 
mains one complication with this proposed experiment. 
YBCO is a bi-layer material, with Y sitting between the 
bi-layers. It is likely that the staggered pattern on the 
bilayers are out of phase, in which case the magnetic 
field at the Y site exactly cancels. A way out of this 
difficulty is to study the 2-4-7 structure where the two 
layers are asymmetric because they are connected to dif- 
ferent charge reservoirs (single chain vs. double chain). 
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It should be possible to have one plane of the bi-layer op- 
timally doped while the other plane (next to the double 
chain) remains underdoped. Obviously, this proposal is 
quite a challenge (but a rewarding one) for the experi- 
mentalist. 

If it is possible to reach H > Hc2, NMR, /i-SR, neu- 
tron, cyclotron resonance and Shubnikov-de Haas experi- 
ments can all be performed to look for the staggered flux 
state. 



V. CONCLUSION 

In summary, the SU(2) formulation of the t-J model 
leads naturally to a picture of the staggered flux phase 
above Hc2 and a stable hc/2e vortex with a staggered 
flux core in the superconducting state. The basic phys- 
ical picture is that the staggered flux state is nearly de- 
generate in energy to the d-wave superconducting state. 
The pseudogap state is described by fluctuations between 
the staggered flux state and the superconductor. It has 
no long range order, but may be characterized by short 
range staggered currents. There may be short range su- 
perconducting fluctuations as well, but these are not de- 
scribed by conventional phase fluctuations alone. As the 
temperature is raised above Tc, the fluctuations initially 
resemble conventional phase fluctuations but gradually 
cross over to fluctuations into the metallic staggered flux 
state, all the while maintaining the energy gap at (0, tt). 
The picture may reconcile the rather conventional x-y 
model behavior observed 30K above Tc with the sur- 
prising persistence of a few vortex-like excitations up to 
1507^. 1^ 

Inside the vortex core, these fluctuations are almost 
frozen out. The core consists mostly the staggered flux 
phase and the tunneling rate between the two kinds of 
vortex can be very small. The small energy difference 
between the staggered flux state and the superconductor 
in the limit of small doping renders this vortex stable. 
This picture suggests a (H, T) phase diagram shown in 
Fig. 4. Below a relatively high temperature scale (of or- 
der Aq ~ ), the spin-gap phase is formed as described 
above. (This was called the staggered-flux phase in the 
SU(2) mean-field theory. |l^) Its onset is a cross-over, 
not a phase transition. Due to the high energy scale, this 
onset is insensitive to magnetic field, consistent with ex- 
perimental findings, Superconductivity onsets below 
a coherence temperature « xt. In a magnetic field, the 



vortex has a core of radius : 



The state inside the core 



forms staggered currents on some slow time scales. At 
Hc2 ~ 2;^, these cores overlap, forming a truly long range 



ordered staggered flux state. This state has a doubled 
unit cell and its Fermi surface consists of small pock- 
ets of area x, consistent with Fermi liquid theory. Thus 
the metallic state generated by a high magnetic field is a 



spin-gap 




FIG. 4. Schematic phase diagram in the {H,T) plane for 
underdoped cuprates. The pseudogap phase onsets below an 
energy scale Aq. This is described by the spin-gap phase 
where the vector in Fig. 2 is disordered. The dashed line is a 
cross-over temperature. The superconducting state appears 
below Tc ~ xt. Its vortex core contains ordered staggered 
currents. For H exceeding Hc2 the vortex cores overlap and 
the staggered flux state is stabilized. 



Fermi liquid state. This state is connected to the pseu- 
dogap phase by an Ising-like phase transition. The long 
range order of the staggered fiux state requires coherence 
among the holes and we expect its transition tempera- 
ture to be xt (i.e. comparable to the superconducting 
Tc) as well. This phase diagram is in contrast to a recent 
proposal by Chakravarty et al. |^^, who suggested that 
the onset of the pseudogap is a genuine transition. In 
their picture the staggered fiux state will extend up to 
the energy scale Ag. The experimental test of staggered 
currents that we proposed should in principle be capable 
of distinguishing their proposal from ours. 

We emphasize that the zero temperature ground state 
in the x-H phase is entirely conventional, consisting of 
a d-wave superconductor, antiferromagnetic insulators, 
and Fermi liquids. At some critical Xc there is a transi- 
tion between the staggered flux state with Fermi pock- 
ets to a Fermi liquid state with a large Fermi surface of 
area (1 — x). The Xc{H) line should terminate at the 
superconducting Hc2{x) boundary. Our picture of the 
zero temperature phase diagram is the same as that pro- 
posed by Chakravarty et al. However, Chakravarty et 
al. asserted that the transition between the two Fermi 
liquids involves a violent change of Fermi surface topol- 
ogy and, by implication, of physical properties such as 
transport measurements. In contrast, we believe that a 
line of continuous transitions with a change of the trans- 
lation symmetry is possible and in fact likely, in view of 
the smooth crossover observed at _ff = above Tc as a 
function of x. The Fermi pockets are elongated and may 
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merge to form a single Fermi surface in the reduced Bril- 
louin zone in the staggered flux phase for x < Xc- The 
restoration of translational symmetry and a large Fermi 
surface can take place continuously by the disappearance 
of the Fermi surface shifted by (tt, tt) which lies outside 
the first reduced Brillouin zone. The scenario of a con- 
tinuous evolution from small to large Fermi surface via 
the "shadow band" was described by one of us some time 
ago.[|0| 

Finally, a third alternative exists; i.e. the staggered 
flux state is never stable (or in other words, it is de- 
stroyed by strong quantum fluctuations even at T = 0). 
In this case, something resembling the spin-gap state be- 
comes the ground state in a high magnetic field and in- 
side the vortex core. If true, this will be the first example 
of a non-Fermi liquid ground state apart from supercon- 
ductivity in dimensions higher than one. Our proposed 
phase diagram offers a very natural route to avoid this 
exotic possibility. 

We would like to stress that even when the vortex core 
is described by the spin-gap state, the hc/2e vortex still 
has a small core energy which vanishes in the x — > 
limit in the mean-field theory. Hence the hc/2e vortex is 
still stable. In fact, the SU{2) vortex is the only mean- 
field theory at present which gives a stable hc/2e vortex 
with a pseudogap in the core. Whether there exists a 
quasi-static staggered current inside the core is a question 
which is difficult to treat theoretically, and which is best 
settled by experiments. 

We end by making a comment on the experiments pro- 
posed by Senthil and Fisher to test for electron frac- 
tionalization. They propose trapping a vortex in a hole in 
a superconductor. When the temperature is raised above 
Tc, the magnetic field escapes, but the Z{2) vortex (vi- 
son) is trapped. Then when the temperature is cooled 
down below Tc, the vison must capture a magnetic flux 
to spontaneously form a hc/2e vortex of either sign. We 
would like to point out that our model of the vortex does 
not exhibit the Senthil-Fishcr effect. While our vortex 
is also a bound state of a magnetic fiux with half a flux 
quantum of the gauge field a^^\ the important difference 
is that the gauge vortex has a finite extent and is not a 
fiux tube. Above Tc, the size of this gauge flux will ex- 
pand to infinity at the same time the size of the magnetic 
vortex does, i.e. the penetration depth of the a'^^^ field 
and the A field both diverge in the normal state. They 
allow the gauge vortex to escape the hole in the normal 
state. 

In principle, the Senthil-Fisher effect, the electron frac- 
tionalization (or the true spin-charge separation), and 
other physics of the Z{2) theory can be readily obtained 
from the SU{2) slave boson theory, if one assume the 
SU{2) gauge symmetry is broken down to Z{2) gauge 
symmetry (which can be achieved by non-collinear SU (2) 
fiux through different plaquettes). |^ With this under- 
standing, the difference between the Z{2) approach and 



our SU{2) approach is clear. In the Z{2) approach, one 
assumes that the SU (2) gauge symmetry is broken down 
to Z{2). While in our SU(2) approach, the SU(2) is only 
broken down to U{1) in the normal state (by a coUinear 
SU{2) flux). The Z{2) and our SU{2) approaches corre- 
spond to different choices of mean-field states of the same 
SU{2) slave boson theory. 
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